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Abstract 



While there does not at this time exist a complete canonical theory of full 
3+1 quantum gravity, there does appear to be a satisfactory canonical quan- 
tization of minisuperspace models. The method requires no 'choice of time 
variable' and preserves the systems' explicit reparametrization invariance. In 
the following study, this canonical formalism is used to derive a path integral 
for quantum minisuperspace models. As it comes from a well-defined canon- 
ical starting point, the measure and contours of integration are specified by 
this construction. The properties of the resulting path integral are analyzed, 
both exactly and in the semiclassical limit. Particular attention is paid to 
the role of the (unbounded) Euclidean action and Euclidean instantons are 
argued to contribute as e"'*^^'/^. 
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I. INTRODUCTION 



Although a complete theory of 3+1 quantum gravity is not yet available, interesting 
arguments for the pair production of black holes the instability of certain vacuum 

states and even for a preferred 'ground state' of the theory f^J^ can be made by analogy 
with field theories that are better understood. A common tool in such arguments is a path 
integral representation of a quantum gravity transition amplitude and/or a related semi- 
classical approximation. Our goal here will be to investigate this representation and the 
h ^ limit by using an analogy with a different class of systems: the 'finite dimensional 
reparametrization invariant models.' 

While they are explicitly defined to have a finite number of degrees of freedom, such 
models possess many of the properties which distinguish Einstein-Hilbert gravity from the 
more common quantum mechanical systems or field theories. They are invariant under 'time 
reparametrizations,' a gauge transformation that mixes coordinates and momenta, and as a 
result they often possess constraints that are second order in momenta. Such 'Hamiltonian 
constraints' are reminiscent of the 'Wheeler-DeWitt equation' of quantum gravity. In 
addition, they share with Einstein-Hilbert gravity the property that neither the Hamiltonian 
function nor the Euclidean action are bounded below. 

Since a bounded Euclidean action is required for common arguments involving analytic 
continuation to Euclidean time, this property has raised concern about how a path integral 
for gravity might be defined and analyzed ||,^. One proposal is to 'rotate the 

contour of integration' in the path integral until the Euclidean action becomes positive 
definite. We will study this question below and compare our results to what would follow 
from contour rotation. 

An important and related point concerns the semiclassical approximation. Recall that 
the Einstein-Hilbert action takes the form 



^~ V~9R + boundary terms (1-1) 

JM 

where g is the determinant of some metric and R is its scalar curvature. One might therefore 
expect that quantum gravity path integrals can be treated as integrals of analytic functions 
and that the contour can be deformed from the original region of integration. Some argu- 
ments for the contribution of Euclidean instantons can be made in this way, by deforming 
Lorentzian metrics to Euclidean ones through the complex plane. However, once deforma- 
tions into the complex plane are allowed, one must in principle include in the semiclassical 
analysis any stationary points which lie within the domain of analyticity of the integrand. 
Because of the branch cut introduced by the ^/g factor in iLlj the metric g effectively takes 
values on a Riemann surface and, at least for appropriate boundary conditions, the station- 
ary points occur in pairs with opposite signs of the action. This point has been made many 
times P^|-p!5| and a similar feature arises in our finite-dimensional reparametrization invari- 



ant models. Which stationary phase points actually contribute to the result will depend on 
the particular contour of integration. This is, however, an issue that must be analyzed and 
one is lead to wonder what determines this contour so that, for example, black hole pair 
creation calculations predict an exponentially suppressed rate for large mass and not one 
that is exponentially enhanced. 
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Our goal here is to investigate these issues in the finite dimensional reparametrization 
invariant context by deriving our path integral from a canonical quantum formalism. The 
existence of a canonical formalism which can be applied to a large class of models and 
which does not require 'deparametrization' or the imposition of gauge fixing conditions 
is a fairly recent development [|T^-|20| whose implications for path integral methods have 
not yet been explored. Such is our aim here, and we will find that the resulting path 
integral has two interesting properties: First, it can be written as an integral over both 
'Lorentzian' and 'Euclidean' paths in which the contribution of Euclidean paths is always 
exponentially suppressed for small h. Second, the Euclidean semiclassical approximation 
yields only exponentially damped contributions; in particular, Euclidean instantons always 
contribute with the weight e"''^-^'/'*, where Se is the Euclidean action of the instanton. 

We begin by describing the aforementioned canonical formalism in section This brief 
review is intended to provide a working understanding of the scheme without addressing all 
of the technical details or supplying all of the motivations. A discussion of these points can 
be found in the literature [|l|,|l3,|19|,|20| . 

After describing the the canonical formalism, we proceed in section ^ to define our 
transition amplitude and to express it as a path integral. This transition amplitude is then 
evaluated in section ^ for two classes of exactly solvable models. In both, it is clear that 
'Euclidean' semiclassical contributions are exponentially suppressed for large \Se\- Section 
V A| rewrites the path integral in a more strongly convergent form which may be useful for 
numerical investigations. Section then [VB| addresses the semi-classical approximation. We 



close with a brief discussion in section VI 



II. THE CANONICAL FORMALISM 



We now describe the specific class of models to be addressed and review the canonical 
formalism of ||l6|JT7|Jl9[| on which our path integral will be based. The goal of this review is to 
provide a working understanding of the scheme and not to describe the technical subtleties in 



detail. For a more complete and rigorous presentation, see |ll6| , p!9| , p0| and for a discussion of 
existence and uniqueness issues, see [^. Essentially the same construction was introduced 
in 



2 1,22 1 in slightly different contexts. 



In this paper, we will study only finite dimensional systems possessing a single constraint 
proportional to the system's Hamiltonian. We take the system to be expressed in canonical 
language in terms of a phase space F, which we shall assume to be T*R" and to be equipped 
with coordinates ,pi,i & 0....n — 1. Finally, we assume that the constraint is of the form 



h = g'^{x)piPj + V{x). 



(2.1) 



As is well known, many such systems arise as 'mini-superspace models' of gravitating systems 
p3|^7| and so form a set of some interest. Note, however, that the method applies to more 
general models [p!6| , p!7| , p!9| , p0| and has been used to construct states of hnearized gravity ||21| 



as well as a Hilbert space of 'diffeomorphism-invariant' states ||T9[ in the loop representation 
approach to quantum gravity. 

The quantization scheme to be followed here is known as the 'refined algebraic method' 
(which is closely related to the 'Rieffel induction method' of |TH]) and may be thought of as 
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an elaboration of the Dirac approach ||2^ , in which the constraints are required to annihilate 



the so-called 'physical states.' Specifically, the refined algebraic approach asks that we first 
quantize the system while completely ignoring the constraint. This provides an 'auxiliary' 
Hilbert space Tiaux in which to work. This space is called auxiliary because it contains much 
more than the 'physical states' that satisfy the constraints. In our case, we will take this 
space to be Haux = L^iH"') with operators X* (coordinates) and Pi (momenta) acting in the 
usual way. For simplicity in our expressions, we define our space using the measure d'^x. 



In this paper we follow the convention of [T^ in denoting classical phase space functions 



by lower case letters while denoting quantum operators by capital letters. We use units in 
which ^ = 1, except in section |V Ij| which explicitly investigates the h ^ limit. 

The next step in the procedure is to 'quantize' the constraint h = 0. For our purposes, 
this simply means that we choose some self adjoint operator H on Tiaux which has the 
function h as its classical limit. The usual ambiguities are present at this level and we 
make no attempt to give a unique prescription. In fact, a somewhat greater ambiguity is 
present here than in quantizing the Hamiltonian of a non-relativistic system. The point is 
that, classically, the constraint h = is equivalent to any constraint of the form f{x)h = 0, 
although such a rescaling can affect our quantum prescription when f{x) is not a constant]^. 
Our viewpoint here is that this is just one more of the many ambiguities that arise when a 
classical system is quantized. 

Now, if the spectrum of H were entirely discrete, the implementation of the Dirac pre- 
scription would be straightforward. Those eigenstates of H with eigenvalue zero would 
become the physical states of our theory and the 'physical Hilbert space' could simply be 
the H = eigenspace of Tiaux- However, in typical cases H will also have a continuous spec- 
trum at zero eigenvalue, for which the corresponding eigenstates will not be normalizable 
in the auxiliary Hilbert space but will instead be 'generalized eigenstates' of H, a kind of 
distribution. 

The strength of the refined algebraic quantization procedure is its ability to form a 
physical Hilbert space from such generalized states using an inner product that is, in a certain 
sense, induced from Tiaux- This has the advantage that any sufficiently 'nice' operator A 
on Tiaux which commutes with H induces a densely defined operator A^^ys on the physical 
Hilbert space 'Hphys- In addition, the map A ^ Ap^yg is a '*-algebra homomorphism,' 
preserving multiplication, addition, and Hermitian conjugation of the operators. Since Tiaux 
is a quantum version of the phase space F, it is through this induction process and the 
auxiliary Hilbert space Tiaux that the *-algebraic properties of the observables on Tiphys are 
connected to the reality properties of the classical phase space functions. This is in fact the 
important point, as the most 'physical' requirement of an inner product is that it give the 



proper adjointness relations to the quantum operators [^,0. This construction is described 



below. The reader is encouraged to consult ||17||r9|j20|| for further details 



We shall in fact assume the spectrum of H to be entirely continuous at if = 0. That this 
case is in some sense sufficient follows from the result |jl9|,^ that the continuous and discrete 
eigenstates of H induce sectors of the physical Hilbert space which are supers elected relative 
to each other. The presence of discrete eigenstates would, however, affect the formulation 



^In particular, the proposals of |2S-^] are not compatible with our choice of inner product. 
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of the path integral in ways that we would prefer to ignore. We therefore content ourselves 
with the observation that many minisuperspace models can be formulated with a constraint 
having only continuous spectrum at if = and restrict to this case; for details see [0 and 
in particular |]l8l for the case of the Bianchi IX model. 



In this situation and under a certain technical assumption concerning the operator H, 
the physical Hilbert space is straightforward to construct. What we would really like is 
to 'project' Tiaux onto the (generalized) states which are zero-eigenvalue eigenvectors of H. 
Of course, since none of these states are normalizable, this will not be a projection in the 
technical sense. Instead, it will correspond to an object which we will call S{H), a Dirac 
delta 'function.' Given the above mentioned assumption on H (see [0), the object S{H) 
can be shown to exist and to be uniquely defined. Technically speaking however, it exists 
not as an operator in the Hilbert space TCaux, but as a map from a dense subspace S of 
'Haux to the (for our purposes, topological) dual S' of S. The space S may typically be 
thought of as a Schwarz space; that is, as the space of smooth rapidly decreasing functions 
on the configuration space. In this case, S' is the usual space of tempered distributions. Not 
surprisingly, this is reminiscent of the study of generalized eigenf unctions through Gel'fand's 
spectral theory |^ and S C TCaux C S' forms a rigged Hilbert triple. 



The key point is then as follows. While generalized eigenstates of H do not lie in Haux, 
they can be related to normalizable states through the action of the 'operator' S{H). That is, 
generalized eigenstates \ipphys) of H with eigenvalue 0, can always be expressed in the form 
S{H)\ipo), where \ipo) is a normalizable state in 5 C Haux- This choice of {ipo) is of course 
not unique and, in fact, we associate with the physical state l^ipphys) the entire equivalence 
class of normalizable states \ip) G S satisfying 

S{Hm = \^p,ys). (2.2) 

Each equivalence classes of normalizable states will form a single state of the physical Hilbert 
space. 

All that is left now is to 'induce' the physical inner product from the auxiliary Hilbert 
space. Naively, the inner product of two physical states \(pphys) and \i^phys) may be written 
{(j)\6{H)6{H)\ip), where |0) and \ip) are normalizable states in the appropriate equivalence 
classes. This inner product is clearly divergent, as it contains [(5(if)]^. The resolution is 
simply to 'renormalize' this inner product by defining the physical inner product to be 

{(l)phys\^phys)phys = (0| 5(i/) |^)a«x, (2.3) 

where the subscripts phys and aux on the brackets indicate the two different inner products. 
Note that ^]3| does not depend on which particular states \ip) G S were chosen to 
represent the physical states \(pphys) and l^pphys)- This construction parallels the case of purely 
discrete spectrum as, if Ph were a projection onto normalizable zero-eigenvalue eigenstates 
of H, we would have [Ph]"^ = Ph- Although S{H) is not strictly speaking an operator, taking 
\(f)) and to lie in S makes the above inner product well defined, as well as Hermitian. In 
the case of the free relativistic particle, this positive definite inner product corresponds to 
the Klein-Gordon inner product on the positive frequency states, but corresponds to minus 
the Klein-Gordon inner product on the negative frequency states. The positive and negative 
frequency subspaces are orthogonal as usual. A similar representation of the inner product 
holds in certain other cases 
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Perhaps the most important feature of this approach is that it first defines the algebra 
of quantum observables (without requiring them to be found exphcitly) and then provides 
a *-representation of this algebra on the physical Hilbert space; that is, a representation on 
Tiphys in which the proper Hermitian conjugation relations hold. From the algebraic point of 
view, this is the fundamental goal of any quantization scheme, and it is this *-representation 
that determines all physical predictions. 

The representation is defined as follows. The observable algebra is defined by the *- 
algebra of observables that commute with the constraint H and act 'nicely' on the dense 
set S (see [0,^] for details). These are the analogues of the (smooth) gauge invariants of 
classical physics. Each such operator A then induces an operator Ap^ys on Hphys through 

Aphysli^phys) = 6{H)A\^), (2.4) 

where again \ip) is any state for which \ipphys) = ^{H)\ip). The operators Ap^ys then satisfy 
the same algebraic and Hermitian conjugation relations as the observables on Tiaux, forming 
the desired *-representation. The use of such operators and the physical inner product 



has been shown to give physically reasonable results in interesting special cases [|T3)0 and 
in the semiclassical limit 



III. A PATH INTEGRAL FOR THE INNER PRODUCT 

Having described how our models will be quantized, we now wish to derive a path 
integral formalism for these systems. To do so, we must first answer the question "Just what 
quantity should we derive a path integral for?" Path integrals are often used to represent 
the 'transition amplitudes' that encode the time evolution of quantum systems. However, 
for the cases we consider, the Hamiltonian explicitly vanishes on the physical Hilbert space. 
Thus, the operator e~*^* is just the identity. 

Nevertheless, we know that the physical states do contain information that we may call 



dynamical (see, for example, [p^ -^ for general comments or |117| , |36[| for a discussion in 
the context of this particular approach). Thus, there should be some mathematical object 
which, more or less, encodes our idea of a 'transition amplitude.' 

Apparently nontrivial path integral expressions for 'transition amplitudes' have in fact 
been studied by a number of authors (e.g., [jl4|,|15|,0). In the minisuperspace context, all of 
these are transition amplitudes between two configurations x and x' (which are the analogues 
of the three-geometries of 3 + 1 gravity), or perhaps their conjugate momenta. As a result, 
we will seek our transition amplitude in the auxiliary Hilbert space Tiaux- 

We take our 'transition amplitudes' to be just the matrix elements of the 'operator' 5{H) 
in Tiaux- That is, we will compute {x\5{H)\x') where |a;) and \x') are generalized eigenstates 
of the coordinate operators X*. Our reasons for this are two-fold. First, while expressed 
in terms of the auxiliary space, such matrix elements contain all of the information about 
'Hphys as they define the physical inner product. Second, when one of the coordinates (say 
x°) is considered to represent a 'clock' and when this clock behaves semi-classically this 
object does in a certain sense describe the amplitude for the 'evolution' of the 'state' x at 
'time' x^ to the 'state' a/ at 'time' x^' . Here x represents the coordinates on a slice through 
the configuration space of constant 
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It is now straightforward to represent this object as a path integraL To do so, consider 
the path integral expression for the operator e~^^^ on Tiaux-, which we expect to exist and 
which can be derived in the usual way by skeletonization (see, for example [^) of paths 
between x at 'time' zero and x' at 'time' N . Note that by 'time' we mean an additional 
parameter that we now introduce; not one of the coordinates x*. We then integrate from 
— oo to oo to turn e~*^^ into S{H). 

The resulting path integral is then 

{x\6{H)\x) = ^ dN J DxDpexp (^i Jj dt[px - h{x{t),p{t))]^ (3.1) 

where dN denotes an integral over the single variable N while DxDp denotes the Liou- 
ville measure on the canonical path space |^3| . Since the configuration variables are specified 



at the endpoints, there is 'one less' set of integrations over the coordinates than over the mo- 

{p\x) 



menta. The notation h{x{t),p{t)) stands for the 'symbol' of the operator 



H, where \p) is an eigenstate of the momenta Pi. 

In the usual way, gauge fixing machinery and redundant integrations can now be intro- 



duced to write 3.1 in a form where its independence of 'gauge' is more explicit. Halliwell 



29| has performed this analysis in the reverse direction by starting with the Fadeev-Popov 



form and introducing the gauge fixing condition = 0. Since he arrives at just our result 



3J] , his work allows us to express the physical inner product in the form 

{x\5{H)\x) = ^J DNDxDp 5(G)A(G')exp [i dt[px ~ N{t)h{x{t),p{t))]^ (3.2) 

where G is now any appropriate gauge fixing function, A(G') is the associated Fadeev-Popov 
determinant, and the sum is over all paths {x{t),p{t), N{t)) in which x'' (t) , pi(t) , N (t) are 
allowed to range over the entire real line. This is the path integral that we shall explore in 
section |V[ For interested readers, the convergence properties of |3]1| and ^]2| are discussed in 
detail in appendix ^ 



IV. EXACTLY SOLVABLE MODELS 

Having derived a path integral for {x\6{H)\x') , the 'physical inner product,' it is of 
interest to see what form this distribution takes in simple cases where an exact analytic 
expression can be obtained. As usual, the cases that we will study are the 'purely quadratic 
ones;' the perturbed Bianchi I model (or free relativistic particle) and the case of coupled 
Harmonic oscillators. 



A. Perturbed Diagonal Bianchi I 

The Bianchi I model is a minisuperspace describing spatially homogeneous spacetimes 
of the form Ai = x H. having a foliation by three-tori with fiat Riemannian metrics (so 
that the tori form spacelike hypersurfaces of A^). In the diagonal version of this model 
the metric is such that, at each point x G M., three mutually orthogonal closed geodesies 
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intersect at x and each encircle an arm of the torus once. This system may be formulated 
on the configuration space Q = with a constraint of the form 



h 



BI 



2 I 2 I 2 
-P0+Pl+P2- 



(4.1) 



In this case, the coordinate x° describes the volume of the three-torus while the coordinates 
and x"^ describe the 'anisotropics,' the ratios of the lengths of minimal curves encircling 
the torus in different directions. We will consider a 'perturbed' and slightly more general 
model on Q = R" for which the quantum constraint is 



-j^ n—l 



(4.2) 



for > 0. Note that without the perturbation (m^), the action 5* would vanish on every 
classical solution and we would be unable to consider the semi classical limit IS"! S> 1. The 
constant factor of | affects none of the results, but conforms to the usual convention for the 
normalization of kinetic terms. 

Such a system looks exactly like the free relativistic particle. However, we will inten- 
tionally avoid referring to this system by that name, as we feel that the physics of the two 
situations is quite distinct. This follows from the fact that the metric which defines the 
constraint's 'kinetic term' has a different interpretation in each of the two cases. A free 
relativistic particle with < is usually interpreted as 'traveling backwards in time,' a 
process that we suppose to be physically disallowed. This leads to the usual preference for 
positive frequency states over negative frequency states. However, in the Bianchi I model, 
a negative po means only that at the x° of the tori is decreasing with proper time - that is, 
that the universe is collapsing. This is not only a physically interesting process but a process 
which classically must occur in some minisuperspace models, such as Bianchi IX Thus, 
we are pleased to include the negative frequency states in our model. 

We now proceed to compute the integral 



{x\6{H)\x') 



1 f°° 

271 J-oo 



dNixle-'^^'^lx') = -5R( / (x|e-*^''|x')) 



(4.3) 



where 3? denotes the real part. The operator e~*^^ is just e"*"^^^/^ times a product of 
propagators for free nonrelativistic particles. As a result, its matrix elements are readily 
seen to be 



(X e 



-iHN\l 



X ) 



ITT /-] n \ 

g 2 ^ 2' 



{2nN) 



-n/2 



exp 



im 



N 



[X — x\ 
im?N 



(4.4) 



for > 0. We can evaluate using 3.471 of [^] to yield, for {x — x') > 



{xm)\x') 



7r(27r)"/2 



{x — x'y 



m 



(l-n/2) 



K(n/2)-i [m^{x - x'Y 



(4.5) 



where K{ni2)-\ is the modified Bessel function (of the second kind) of order (n/2) — 1. 
Similarly, for (x — x)^ < 0, we find 
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{xm)\x') 



(27r)"/2 



(X — X' 



'\2 



m 



(l-n/2) 



COS 



n — 1 



J(n/2)-l [mJ-{x - X'Y 



sm 



vr 



n — 1 



iV(n/2)-i [m^J -{x - x'y 



(4.6) 



An intuitive feel for just why the answer is of this form may be obtained by noting that it 
describes a certain correlation function in free scalar field theory, or alternately by performing 
this integral in the semiclassical approximation. For the interested reader, the semiclassical 
analysis is carried out in appendix 0. 

^2 



Note that, for 



[X — X 



^ m^, the matrix elements are roughly sin(m 



[X — X 



/^2 



) or 



cos mw— X — X' 



'\2\ 



depending on the number of degrees of freedom. That is, they include 
equal contributions from what might be called the 'positive and negative frequency parts.' 
Nevertheless, when 
nential exp(— m 



[X 



X 



i\2 



X — X 



i\2\ 



^ m , the matrix elements contain only the decreasing expo- 
This occurs even though the Euclidean action is unbounded 
below. Similar results hold for the case im? = —k"^ < for which the resulting physical 
inner product is obtained by replacing m with k and 



X 



X 



l\2 



with 



\x 



X 



'Y in 0| and 



As might be expected, this is related to which points of stationary phase contribute to 
{x\5{x)\x') in the semiclassical approximation. We will return to this point in section and 
appendix ^. 

It is interesting to compare the results [4.5| and [4.6| to what one might obtain by a 'contour 
rotation prescription' of the sort suggested in The idea is to rotate the contour of the 
conformal mode [represented here by x^] to imaginary values [x° — > ix^] so that the Euclidean 



action becomes bounded below. We can certainly perform the integral [4.3| taking —p^ to be 
replaced by pg in the Hamiltonian h. The result is just [4.5| above with (x — x'Y replaced by 
{x — x')\ = I]"=o(^*)^' since (x — x'Y^ > 0. The task is then to analytically continue back; 

x')^ > 0, this clearly yields |4.5|; the correct answer from our point of 



X 



-tx 



For 



X 



view. However, for 



X 



X 



/\2 



< 0, the analytic continuation is ambiguous. Because of the 



branch cut that defines the square root, the answer will depend on which path is followed 
through the complex plane. As can be seen from [4.6| , what we would call the 'right' answer 
results from combining these two contributions with equal weights and with a phase that 
depends on n. 



B. Coupled Harmonic Oscillators 

In order to more thoroughly explore the properties of exactly solvable models, we now 
study the case of coupled harmonic oscillators. However, to fit with the general character 
of minisuperspace models, we will take one of the oscillators to have negative energy. That 
is, we will again take our system to be defined on the phase space F = T*R", with the 
constraint 

o = h = h-pi - ix^co^ + Y.{pi + i^r-^')] ^kp' + -^'A (4.7) 

i=l ^ 

for uj'^ > 0. We could, in principle, also study the case where one or more of the Harmonic 
oscillators is inverted (with Hamiltonian |(p^ — tu^x^)), but we will not do so here. For 
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simplicity, we assume n to be even. Models of the form [4.7| do arise as minisuperspace 



models of gravity interacting with scalar fields |^ and many of their aspects have been 
studied in the literature |4^-^ . 

Strictly speaking, however, this model does not fall into the class allowed by section 
fl[ Since it is a sum of harmonic oscillator Hamiltonians, instead of the constraint H 
having purely continuous spectrum at zero (when quantized in the usual way on L^(R"')), 
its spectrum is purely discrete - not just at zero, but everywhere. This creates a number 
of subtleties, such as the fact that there are rather few zero eigenvalue eigenstates of the 
constraint unless all of the oscillators have commensurate frequencies. Furthermore, unless 
an appropriate constant is added to iJ, it will in general have no eigenstates with eigenvalue 
zero. Thus, a better way to quantize this model might be to rewrite the constraint (see, for 
example, |T7|JT8| , |36[] ) in such a way that the spectrum is purely continuous^]. Unfortunately, 
this necessarily destroys the purely quadratic form of the action and thus the exact solvability 
of the model. 

Our purpose in studying this model is not to examine the detailed predictions of the 
cosmological scenario of EH], but rather to investigate the general mathematical properties 



of the expressions and p.2| for the physical inner product. Thus, we will quantize this 
model using the constraint operator 



n-l 



H=-\ 



-In 



(4.8) 



The frequencies and constants in [4.8| are specifically chosen so that there seem to be a suffi- 
cient number of zero-eigenvalue eigenvectors. We shall not concern ourselves with whether 
or not such choices are 'natural.' 

However, if we are to proceed in this way, expressions and p.2| must be modified. This 



is because their aim was to calculate the operator biH^ and, now that the spectrum of H is 
discrete, this object is highly divergent: when acting on a normalizable state I?/') for which 
= 0, S{H) cannot possibly be defined. 
As stated in section the analogous object for the case of discrete spectrum is the 
projection Ph=o onto the zero eigenvalue subspace for the operator H. An expression for 
{x\Ph=o\x') analogous to and ^]2| can be found by realizing that the 'evolution operator' 
^-iNH jj-^ ^Yii^ case, periodic in with period 2tt/uj. This of course is the source of 
the would-be divergence in 5{H) = ^ J^^dN e~*^^, but it also allows us to express the 
projection Ph=o as the integral over a single cycle of A^: 



/vr J-n/w 



■iNH 



UJ 



TC Jo 



n/uj 



dN ( 



e + e^^^ 



(4.9) 



Similarly, the analogues of |3.1| and |3.2| are given by taking A^ to live on a periodic interval 
of length 271 /u. 

As in the perturbed Bianchi I model, the operator e~*^^ is a product of evolution 
operators for nonrelativistic particles which can be evaluated exactly; perhaps most easily 



^This continuous spectrum form of the constraint is in fact in better accord with the global 



properties of the original gravitational model of |46| 
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by using the semiclassical approximation (which is exact for quadratic Hamiltonians) . For 
- > > 0, the result is 

(x|e-*^^|x') = ( ^^y^'e^f{i-5)e*5=K-y;^) (4.10) 



27r sincuA^ 



where 



n r /An UJ {x^ + X ^) COSuN - 2x ■ x' 

Sci(x,x;N) = —- (4.11) 

^ ' ' ^ 2 sincjA^ ^ ^ 

is the action of the least action path between x and x' traversed in the given time N. Note 
that this path is unique for N 0, tt/u. Here, x ■ x' denotes —x^x^' + YaZI x'x'' . 



All that remains is to evaluate the integral over A^ in [4.9| . For the case n = 2, we may 
refer to More generally, if we define a = —uj{x'^ + x'^)/2 and h = ux ■ x' then the 

result may be expressed as 

\ (n/2)-l 

{x\Ph=o\x') = 2cu3fJ ^ (27r)-("/2),gif (i-t) L ] H^^\Vlf^) \ , (4.12) 

r{i) 



using 6.677 from [Q. Here, Hq is the usual Hankel function of the first kind and the square 



root is defined by \fre^ = r^/^e*^/^ for < ^ < vr. Thus, when > b"^ 4.12 takes the form 



{x\Ph=o\x') = — (27r)-"/2 ^ K,iV^^^). (4.13) 



TT 



The results [4.12| and [4.13| may in turn be written as a sum of Hankel functions of orders 



0<iy<n/2 — Ihy using the usual Bessel function recurrence relations. 



The matrix elements |4.12| have the same general structure as those of and |4.6| . That 
is, they are a Bessel function of order n/2 — 1 of the square root of a function of x and 
x'. Recall that in the case of perturbed Bianchi I, the behavior of this Bessel function 
was determined by the sign of {x — x')^; that is, by whether the points x and x' could be 
connected by a 'Lorentzian' classical solution (with real lapse A^) or by a 'Euclidean solution' 
with imaginary lapse. A similar phenomenon occurs for our coupled oscillator model. Pairs 
of configurations {x, x') may be separated into three distinct classes: those for which — | > 1, 
those for which 1 > — | > —1, and those for which —1 > Only for the first category 
(— I > 1) does a 'Euclidean' classical solution (with imaginary lapse) exist. Such solutions 
come in pairs with Euclidean actions Se = ±\/a^ — b"^- Similarly, a 'Lorentzian' solution 
(with real lapse) exists only for the second case, 1 > — | > —1. These solutions also come 
in pairs, with actions S = rLy/b"^ — a?. For the third case (— | < —1), the configurations x 
and x' may be connected by complex classical solutions with 3? A^ = - (but no such solutions 
exist for the first two classes of pairs). The Euclidean action is also real for this last case, 
and is again given by Se = ±V«.^ — 6^. 

The semiclassical approximation is valid (in our units) when 15*1 = \J\b'^ — a?\ ^ 1 so 
that, once again, when x and x' are connected by a Lorentzian solution, this approximation 
contains contributions both of the form e*'^ and e~^^ which are equally weighted up to a 
phase. However, when the connecting solution is Euclidean (or complex with real Euclidean 



11 



action), only the exponentially decreasing solution contributes and the leading semiclassical 
term is e"'"^^'. We follow the usual convention of iS = —Se so that this exponentially 
decreasing factor corresponds to the stationary point with positive Euclidean action. 

This seems to indicate a general property of the semiclassical approximation (and perhaps 
of the entire expression) for and |3.2| which we shall investigate further in section |V[ Note 
that this case is more subtle than that of the perturbed Bianchi I model since that model 
had the property that (for > 0) the Euclidean action of a Euclidean solution is always 
positive when the 'Euclidean lapse' iN is positive. In contrast, for the coupled harmonic 
oscillators, a Euclidean solution with positive Euclidean lapse can have a Euclidean action 
with either sign. This, therefore, is much closer to the generic case. 

As in the perturbed Bianchi I model, we briefly compare [4.12| with a contour rotation 
scheme where is rotated to ix^. Again, this gives just the right result in the 'Euclidean 
sector' (a^ > 6^), but the branch cut of the square root creates an ambiguity when ana- 
lytically continuing back to the 'Lorentzian sector' < b"^. Our result 4.12 includes both 
possible outcomes, equally weighted up to an n-dependent phase. 



V. THE GENERAL CASE 



In sections [IV A| and [IV B| , we discoverd several interesting properties of the exact re- 
sults 



a and 4.12 



We will now argue that these properties should hold in general. 
Specifically, we first show in section [V A[ that the inner product may be expressed as a path 
integral over both 'Lorentzian' and 'Euclidean' paths (and combinations thereof) in which 
the contribution of the 'Euclidean' paths is exponentially suppressed for small h. Such an 
expression is more convergent than ^]T| and and may be of use for numerical investiga- 
tions. A representation of this form also indicates that Lorentzian instantons (with either 
sign of the action) and Euclidean instantons with positive Euclidean action will contribute 



to our transition amplitude in the semiclassical limit. We then argue in section [VB| that, in 
addition, Euclidean instantons with negative action do not contribute to the semiclassical 
limit. 



A. A simplified expression 

Our general strategy will be to simplify the path integral by performing the integrals in 
3^ over the momenta and over the lapse; thus, we assume that we may change the order 



of integration. The result will (almost) be a configuration space path integral in which, in 
effect, the constraint equation 

9^A^)^ + Vi^) = (5-1) 

has been solved for the lapse N. That is, we will obtain a path integral based on a 'Baierlein- 
Sharp- Wheeler-like' form 



5* 



J dt ^JgijX^x^V{x) (5.2) 
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of the classical action. 

We begin with the expression |3.2| for the matrix elements {x\S{H)\x') . Having written 
the integral in a form where Fadeev-Popov technology applies, we are now free to use any 
gauge fixing condition we wish. Because we would like to integrate out the momenta, we 
will use a gauge condition of the form G{x) = where G is a function of the coordinate 
variables only (see for a discussion of canonical gauges in this context). In addition, we 
assume VG ■ VG > 0; that is, that VG is 'spacelike.' Unfortunately, G = is not really 



a 'good' gauge condition ||5^ as there will always exist points in phase space at which the 
Poisson bracket {G, H} = p-VG of G with H vanishes (such as where all momenta vanish). 
Nevertheless, such a condition is often 'good' on all but a set of measure zero. Having 
admitted our treatment to be heuristic, we assume that we may use such a gauge condition 
below. 

Now, since there is only a single constraint, the Fadeev-Popov determinant takes the 
particularly simple form of a product over times t: 



A{{G,H})=1[ I {G, H}it) I = n I VG(t) Pit) I 



(5.3) 



where the inner product -t denotes a contraction through the (co) metric g^^{x{t)). 

For common factor orderings of the quantum constraint H, the symbol h{x,p) will be 
quadratic in the momenta, p. We shall assume that it takes the form 



h{x,p) = g'^{x)piPj + Vg{x) 



(5.4) 



where the subscript q on V indicates that the potential may receive 'quantum corrections' 
such as terms proportional to the curvature of gij []5l| , p5| . In order to be explicit, we shall 
assume the signature of g to be everywhere (1, n — 1). 

Most of the the momentum integrals are of the form dp e^iv^-^P ) and can be per- 
formed exactly using stationary phase methods. However, the presence of the absolute value 
I VG ■ p\ means that, for each value of t, there will be one momentum integral which cannot 
be performed in this way. It will be sufficient for our purposes to leave this integral un- 
done, and to simply perform the others. The component of p in the direction of VG(a;(t)) 
will be denoted p\\. Similarly, n\\{x{t)) will be the unit vector in the direction of VG(a;(t)). 
Performing these momentum integrals yields the expression: 



{x\5{H)\x') = J DNDxDp\\exp 



I- I dtgi 



x^x^ 



NiV,ix)+pl 



, 7rsign(iV) / n-3 \ 1-n 

t 



detg^ 5[G{x{t)] |VG(x(t))| 



(5.5) 



where the usual product of normalization factors at each time has been dropped (they may 
absorbed into the definition of the measure). Here g:^ is the induced metric on surfaces of 
constant G(x). 

Finally, we wish to perform the integrals over A^. As in the explicitly soluble models 
discussed in section these integrals may be expressed in terms of Bessel functions. Note 
that each such integral is of the form 



. 7rBign(JV) / n-3 

dNe 5 I — 



^)|iV(t)|(l-")/2 



exp 



1 x^x^ 



N{V,{x)+pI 



(5.6) 
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This is essentially the same integral with which we were faced in section |IV A| . Introducing 
V = V + n'^^Pp k = gijx'^x\ and dropping constant normalization factors, the integral yields 



when kv > and 



-A:/.)(^-)/^e-(^)<_3)/,(v^) 



(5.7) 



(5.8) 



when kv < 0. 

As a result, the physical inner product may be written 



{x\S{H)\x') = f DxDpiiYl Iik,v) 6iG)\VG\\p\\\^-detg- 



(5.9) 



The connection to the Baierlein-Sharp-Wheeler-like ^.2| should be clear: since Bessel func- 
tions are exponentials at large arguments, if the momenta p\\ were replaced by their semiclas- 



sical values x ■ VGyv/k then |5.9| would be a path integral based on the action ^.2|, although 
with a rather complicated measure. 

This expression may be thought of as a sum over both 'Lorentzian' and 'Euclidean' bits 
of path at each time. Here, we say that a segment of path is Lorentzian when kv > 0, so 
that a real lapse would be associated to this segment by solving the constraint 



_k_ 

iV2 



+ v = 0. 



(5.10) 



Similarly, a bit of path with A;/f < is 'Euclidean' as the corresponding lapse is imaginary. 
Note that this occurs despite the original expression 3J. being an integral only over real 
classical lapse A^. As in our exactly solvable models, the Lorentzian bits contribute with 
both signs of the action but the Euclidean path bits are always exponentially suppressed when 
the corresponding \Se\ is large. A similar result is obtained if other momentum integrals 
are left undone as well, so long as the momenta at each time are integrated over at least one 
two-plane with signature (—,+). 

The argument of every Bessel function is now manifestly positive. As a result, the fact 
that neither the Lorentzian nor the Euclidean action is positive definite has disappeared 
from sight. If we define a slightly modified function by Ie{k, v) = I{k, v) for kv > and 



-k/vr-y^e-("-^)HlZ,y,i[l - ^\e\\^v) 



(5.11) 



for kv < 0, then I^{k,v) vanishes in the limit of large kv (for e ^ 0). Our matrix elements 
may therefore be expressed in the form 



(5.12) 



{x\6{H)\x') = lim / DxDpnYmHt),v{t)) 

X 5{G{x{t)) \VG{x{t))\ bill ^-deig^ 

for which the integrand is exponentially decreasing. As it converges more strongly than ^TT] 
we may hope that |5.12| will be of use in numerical computations. 
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B. Instantons and the Semiclassical Approximation 



The expression |5.9| for our transition amplitude includes an explicit sum over both 
Lorentzian and Euclidean paths, as well as arbitrary combinations of the two. In the semi- 
classical approximation, this result indicates that both Lorentzian and Euclidean stationary 
points contribute to physical inner product, so long as the latter have positive Euchdean 
action. In fact, since the contribution of Euclidean path-bits is exponentially suppressed 
in |5.9| for small h, it appears that stationary points with negative Euclidean action never 
contribute, as was suggested by the examples of section |V[ However, we have not yet 
shown this carefully. In particular, it is possible that a stationary point of the action, while 
not actually on the contour over which the integral is performed, might contribute to the 
semiclassical approximation^. In general, in fact, we do expect such stationary points to 
contribute, as the contour of integration could easily be deformed to reach a complex sta- 
tionary point lying just off the real axis. If the contour, say for the variable k in |5.9| , can be 
deformed far enough, it could wrap around the branch point at A; = to reach a stationary 
point on the 'second sheet' of the Riemann surface on which k lives. Such a stationary point 
would then contribute as e"'"''^^' to the inner product {x\S{H)\x'). We will now argue that 
this does not occur. Because we study the semiclassical approximation, the factors of h will 
be restored below. 

Let us first note that, if the semiclassical approximation is to hold, our matrix elements 
must be of the form 

{x\6{H)\x') =e'^("'"')/^C(x,x') (5.13) 

where C{x,x') is slowly varying in comparison with e'^six,x')/h^ particular, since S{x,x') 
is continuous on regions TZ of R" x R" where the number n of connecting paths in constant, 
C{x,x') must be continuous there as well. 

In general, we expect the union of the boundaries of such regions TZ to have measure 
zero so that the union U of the interiors has measure one. We therefore assume that this is 
so and consider only the open subset U on which {x\6{H)\x') is continuous. 

If the matrix element {xo\S(H)\xq) takes the form e+l'^^^^'^'^o)!/^ (to leading semiclassical 
order^) at some pair {xo,Xq) G U, it must in fact take this form on some open rectangular 
region Vn, x Wn (Vn C R", Wn C R") containing (xq, x'q). If the diameters of Vn, Wn are much 



less than h 



-1 



and h 



dx 

VhXWn- As a result, {x\6{x 
where (say) A = ^\Se{x,x')\. 



dSE 



dx' 



respectively, then S{x,x')/h is essentially constant on 



x') is essentially real on V/j x Wn, and satisfies {x\6{H) \x') > e 



We have already seen an example of such behavior as the Euclidean paths that are explicitly 



included in 5.12 correspond to stationary points of the N integrals at imaginary lapse (and therefore 



off the original real contour). 

^The fact that a Euclidean path dominates does not necessarily require the matrix elements to 
be real and positive; there may in fact be an overall phase that we have neglected. The point is 
that this phase is both independent of H (for small h) and slowly varying with x and x'. 
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We now choose two states |0) and in 5 whose representations (a;|0) = f{^^^), 
{x\ip) = g{^^^) are positive real, supported in Vn and Wn respectively, and have /(?/) and 
g{y) independent of h. Let us define 

= ^""/Rn(x|0)d"a; = /r„ f{y)d''y, 
H = ^""/an {x\ip)d'^x = Jnngiy)d''y 
h = ^""/r" = |(a;|^)|2d"a;. (5.14) 

It follows that a^, a^, b^j,, are real, positive, and independent of h. Note that the physical 
inner product of our two states satisfies 

{<l)n\S{H)\^n)>fi^''a^a^e^/^ (5.15) 

to leading semiclassical order and that the expression on the right diverges a.s h —>■ 0. 

We will now derive a contradiction. To do so, we return to our original expression for 
the inner product 

r+oo 

(0,|5(iJ)|V^,)= / (0,|e-^^^/^|^,). (5.16) 



Recall that, for {(pn), lipn) G S, this integral converges at large N. Effectively, this is because 
the states \(pn), \i^n) are characterized (to accuracy e) by some minimal energy scaleQ En 
and the integration over the region |A^| > ^ h/E^ yields a negligible contribution. Note, 



however, that since \4>h)i iV'ft) have norms yh^b^ and yh^b-^ and since the operator e ^^^/^ 
is unitary, we have the bounds 



{Me-'^'^'^^m <K''^b^b^ (5.17) 

and 

{MS{H)\i:n) < h^'^'C^/b^/En (5.18) 

where C is some constant independent of h. In order to compare this bound with the 
semiclassical expression |5.15| , we estimate E^ as follows. 

Recall that h = g''^{x)piPj + V{x). The states {(pn) and {ipn) are characterized by coordi- 
nates Xq and Xq (which are independent of h) and by momentum scales 



h 

Pn 



diam(V7j) 



dx 



h 

Ph 



dx' 



(5.19) 



diam(iy;i) 

which are also independent of Ti. As a result, E^ ^ g{xo)Ph + ^(^^o) is independent of h 
and the bound |5.18| vanish /i — > 0, in contradiction with p.l5| . Here, g{xo) is some 

bound on the components of the g^^{x) at x = Xq. We thus conclude that, at least for 
minisuperspace models, Euclidean instantons always contribute to our transition amplitude 
as e~''^^'/^. Similarly, complex instantons must contribute as e*^'^/^e~'^'^/^', where 3 denotes 
the imaginary part. 



-"This is the energy at which 



{<l)n\5{H)\^n) - mSiH - En^n) 



< e. As a result, it is related 



to a sort of continuity |2C] of the spectral representations of {(pn) and \iph)- depends on the 
functions / and g as well as on h and the accuracy e. In particular, it is not a property of the 
Hamiltonian H alone and in no way corresponds to a mass gap for the system. 
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VI. DISCUSSION 



The goal of this work was to use a canonical formalism (the refined algebraic or Rieffel 
induction method) to represent a minisuperspace transition amplitude as a path integral and 
to investigate this expression both exactly and in the semiclassical limit. We have seen that 
the path integral may be written in the form |5 .121 which explicitly sums over both 'Lorentzian 
path-bits' (which appear with both signs of the action) and 'Euclidean path-bits' which are 
exponentially suppressed. As a result, we conclude that any Lorentzian or positive action 
Euclidean stationary point may contribute to our transition amplitude in the semiclassical 
limit. We also argued in section |V B| that Euclidean instantons with negative Euclidean 
action (which would be exponentially enhanced) do not contribute in the semiclassical limit. 

It seems reasonable to assume that similar results hold in full 3-1-1 gravity, Kaluza- 
Klein theory, dilaton gravity, or any other diffeomorphism invariant theory of gravity with 
indefinite Euclidean action. However, it should be noted that the most interesting instantons 
in quantum gravity (e.g., [p]-|6|,p!3|- p^p6| -|58|] ) are rather far from our minisuperspace models. 
In particular, they involve processes in which the spatial topology of the initial state differs 
from the spatial topology of the final state; an effect which appears to be ruled out by 
construction in our context (see however, [|l^,^). An argument could certainly be made 



that such instantons are qualitatively different and that our results have no bearing on their 
interpretation. On the other hand, we note that the arguments of [VB| did not rely on 
the details of our model, but followed from the unitarity of e~*^^ and the existence of an 
appropriate subspace S C Haux- One would expect such arguments to generalize readily 
if refined algebraic/Rieffel induction methods are applicable at all. At the very least, our 
arguments are suggestive, and it is worthwhile to briefly discuss below the hypothesis that 
our results do generalize to such cases. 

Many instanton calculations, such as the pair creation calculations of [0-0], implicitly 
assume that the relevant stationary point is the one with positive Euclidean lapse. So 
long as the corresponding Euclidean action is positive (as it is in all of the pair creation 
examples [jl|-0]), this is just the conclusion that would follow from section |VB . We will 
therefore concentrate on situations where the Euclidean action is negative. Perhaps the 
most interesting use of such instantons was seen in the arguments of Baum Hawking 
and Coleman that the cosmological constant A should vanish. They supposed that, 
in some way, the quantum state of the universe provides a probability distribution for A and 
proceeded to estimate this distribution through instanton calculations. The large negative 
action of four-sphere instantons for small positive A was used to argue that this distribution 
is exponentially large (or even larger [^) for A near zero. 

While we have been interested in the transition amplitude and not in a particular quan- 
tum state, we note that every physical state {ipphys) in our scheme may be expressed in the 
form 



x\i>phys) = / dx'{x\5{H)\x'){x\ 



(6.1) 



for some G 5. As a result, the arguments of section |VB| imply that no physical state has 
a 'wavef unction' {x\^phys) which behaves as e''^^*^^^' in any region of superspace. This sug- 



gests that the instantons of |56-58| would not in fact contribute to the desired distribution. 
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Following a similar line of reasoning, arrives at this same conclusion. Other arguments 



against and comments on the Baum-Coleman-Hawking mechanism include [|60H 64 

Another point to be addressed is the argument of |13[ that only stationary points with 
"^■s/Qe > (where is the Euclidean metric) should contribute to the path integral. In our 
notation this restriction would be ^{iN) > 0; i.e., positive Euclidean lapse. The argument 
of was based on considering the normalizability of an induced quantum state for the 
matter (non-gravitational) fields. The issue concerned whether this state is an exponentially 
growing or decaying function of these fields; that is, it had to do with the state functional 
at large values of the matter fields. While we have argued that stationary points with 
3fJ(iiV) < may contribute, the fact that the Euclidean action of matter fields is positive for 
^{iN) > means that for large values of the matter fields our condition of positive Euclidean 
action is equivalent to the condition 3?(iiV) > of |13[. This is as one would expect since 
exponentially increasing functions do not define generalized eigenstates of the constraint. 



As a result, we see that our arguments are consistent with the general requirements of ||T3 
for reproducing quantum field theory in curved spacetime. 



ACKNOWLEDGMENTS 

This project sprouted from a conversation with Andrew Chamblin, to whom the author 
would like to express his gratitude. The author would also like to thank Andrei Barvinsky, 
Jim Hartle, Gary Horowitz, Joe Polchinski, and especially Jorma Louko for helpful comments 
and discussions. This work was supported by NSF grant PHY95-97965. 



APPENDIX A: ON THE CONVERGENCE OF THE PATH INTEGRAL 

REPRESENTATIONS 



In this appendix we discuss in detail the convergence properties of the path integrals |0 



and ^]2|for the physical inner product. These expressions do in fact converge, so long as all 
integrations are appropriately interpreted. 

Let us begin by noting that {(ple'^^^lip) is well-defined for all \ip) E Haux, and 
in particular for all {ip) G S. This follows from the fact that e*^^ is a unitary (and 
therefore bounded) operator on Haux- As a result, {x\e~^^'^\x') is a well-defined distribution 
on iS X iS; that is, (x|e~*^^|x') is a member of the dual space S' x S'. 

Note that this holds for all real A^, even though the expression for {x\e~^^^\x') (for 
the case H = —Pq + = m^, for example) as a function of A^, x, and x' may have an 
essential singularity at A^ = (see, e.g. ||l^JT5[| ). This is nothing more than the fact that 
the distribution 6{x — x') cannot be represented by a smooth function. In fact, {x\e~'^^^\x') 
is a continuous function of A^ in the topology of S' x S', even at A^ = 0. There is thus no 
difficulty with the A^ integral at any finite value of A^. 

We expect the result of this integral to be {x\S{H)\x') . Our technical assumptions guar- 
antee that the limit 

lim f dN{x\e-'"^\x') (Al) 
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does in fact converge to {x\S{H)\x') (in the topology of S' x S'), and does so independently 
of how the limits of a and b are taken [^]. Any lack of convergence of |3.1| must therefore 
arise from the path integral representation of (x|e~*^^|x'). It is to this expression that we 
now turn. 

Our assumptions (from |2^) guarantee that (l — iHt) maps S into S so that the operator 
{1 — iHN/kY is also a member of It follows that the integrations in the k-skeletonized 

path integral 



must converge in the sense of defining elements of S' x S'. 
In addition, for fixed states Ich) G S and \ih) G 5, we have 



qm] (A2) 



{<j>\e-n^) - {<j>\{l - zHtm < Ct' (A3) 

for small t and some (state-dependent) constant C. Thus, the k ^ oo limit of the k- 
skeletonized path integral converges (as a sequence in S' x S') to (x|e~*-'^^|x'). We conclude 



that, when appropriately interpreted, the expression 3.1 does in fact converge, despite the 



oscillatory nature of the exponential and the fact that the 'Euclidean action' is unbounded 



below. Since 3.2 is designed to exactly reproduce 3.1, we have established the convergence 



of this expression as well. 



APPENDIX B: STATIONARY POINTS AND ANALYTICITY OF THE 

MEASURE 



We now show how the the analytic structure of the measure in |4.3| determines the overall 
form of the results [1.5| and |4.6| by evaluating this result in the semiclassical approximation. 
The structure of the results [4.12| and |5.9| is determined in essentially the same way. 



Consider the integral 4.3: 



mH)\x') 



— / dN 
1 

2^ 



gif(l-f)sign(iV) 



(27r|Ar|)V2 

g— i7m/4 



exp I 



A^ 



{X — X 
711? N 



l\2 



dN , 

(27rA^)"/2 



g(j7r/2)sign(Af) 



,m 



l\2 



N - 



{X — x\ 
imP'N 



where the function N^/'^ is (-\/iV)" and \/N is defined to have a branch cut along the positive 
imaginary axis. Note that sign(A^) may be interpreted as an analytic function which is 
constant on the left and right half planes with a cut along the entire imaginary axis. 

Now, for large A^ (and im? > 0), the exponential factor decays in the lower half plane 
but grows in the upper half plane. As a result, the contour may be closed in the lower half 
plane as shown below. 
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o 



o 



A 



\1/ 



± 



when 



/\2 



On this diagram, the filled black circle represents the essential singularity at 
open circles are the Euclidean stationary points at 

the *'s are the Lorentzian stationary points at = ± 



-{x—x')^ 



= 0, the 
> 0, and 



-{x—x')'-^ 



for [x—x') < 0. The vertical 
dotted line is the cut that defines sign(A^), while the horizontal line represents the real axis 
- the original contour of integration. The dashed line going up and then back down the 
negative imaginary axis is a new contour to which the original one may be deformed. Note 
that, because of the cut, the integrations along the negative imaginary axis do not cancel, 
but instead give equal contributions and add together. Since, for {x — x'Y > 0, the integrand 
vanishes as approaches zero from the lower half plane, the essential singularity at A^ = 
does not contribute in the semiclassical limit. This limit is therefore dominated by the 



stationary point in the lower half-plane, for which the integrand in eq. ^1] is exponentially 
suppressed by a factor of e in agreement with the leading behavior of 
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